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Abstract 

Within the context of supersymmetric theories, explaining a 125 GeV Higgs motivates a con- 
sideration of a broader range of models. We consider a simple addition to the MSSM of a "Sister 
Higgs" (Sd), a Higgs field that participates in electroweak symmetry breaking but does not give 
any direct masses to Standard Model matter fields. While a relatively minor addition, the phe- 
nomenological implications can be important. Such a field can be naturally charged under an 
additional symmetry group Gg- If gauged, the Higgs mass is naturally much larger than in the 
MSSM through an NMSSM-type interaction, but with playing the role of H^- The addition of 
the sister Higgs allows new R-parity violating operators HaH^E, which are less constrained than 
conventional leptonic R-parity violation. Considerations of unification motivates the presence of 
colored G^-charged fields. Production of these G-quarks can lead to new b-rich final states and 
modifications to decays of gluinos, as well as new opportunities for R-parity violation. Unlike a 
conventional fourth generation, G-quarks dominantly decay into a light jet and a scalar (poten- 
tially the Higgs), which then generally decays to b-jets. The presence of additional sister charges 
allows the possibilities that lightest sister-charged particle (LSiP) could be stable. We consider the 
possibility of an LSiP dark matter candidate and find it is generally very constrained. 
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I. INTRODUCTION 



As data continue to arrive from the LHC, the simplest models of supersymmetry have 
become increasingly constrained. In particular, the recent discovery of a Higgs-like feature 
at 125 GeV [U E] is tantalizingly close to the tree level prediction of the MSSM, but is large 
enough to require large radiative corrections, and pushes the MSSM into narrow corners of 
parameter space [3J. At the same time, no clear evidence for any of the other superpartners 
has appeared. Both of these results can be accommodated by a higher SUSY scale, but at 
the cost of increasing the tuning of electroweak symmetry breaking scale. 

Alternatives for raising the Higgs mass have been varied: the simplest possibility maybe 
the NMSSM, in which a singlet S is added to the MSSM together with the superpotential 
term XSHuHd- The resulting quartic coupling \huhd\^ raises the Higgs mass at small tan/3. 
However, aside from the aesthetic issues with including a complete singlet, achieving a 
large enough quartic requires a superpotential coupling which must be near the limit of 
perturbativity at the GUT scale. At the same time, if a much larger value for the quartic 
were possible, then improved possibilities for naturalness arise |3j. 

Large quartics can be natural, for instance if the singlet 5* and/or Higgs is a composite 
[U |5] (alternatively the new states can be integrated out yielding quartics from Kahler 
potential terms [6]). New matter at high energies can boost the quartic somewhat by raising 
the standard model gauge couplings in the UV [H [S]. 

While the D-term quartics are fixed by supersymmetry, they can be enhanced by non- 
decoupling D-terms from an additional gauge group P, [10]. While appealing in principle, 
most models require significant additional matter content to maintain unification. 

Extending the MSSM to include new fields is dangerous however, as they can introduce 
new sources of FCNCs or significant corrections to precision electroweak observables. There 
is a simple addition, however, that can evade these concerns: namely, the inclusion of 
additional Higgs multiplets, which we denote S„ and S^. These fields carry the same SM 
quantum numbers as Hu and Hd-, respectively. To avoid FCNCs these fields should have no 
tree-level, renormalizeable couplings to SM fermions. At the same time, there is no reason 
that they cannot participate in EWSB by acquiring significant vevs. Such fields we refer to 
as "sister" Higgs fields. 

The inclusion of additional Higgs fields is certainly not new. What we are describing is 
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largely the generalization of the Type I 2HDM to a supersymmetric context. Similar non- 
coupling Higgses were employed to construct R-symmetric SUSY theories [TT]. As we shall 
explore in this paper, however, there are important phenomenological consequences when 
allowing fields such as these in the low energy theory. 

First, an important point is that for NMSSM-like enhancements of the Higgs mass, the 
coupling SHyT^d is just as effective as SHuHd, i.e., the object that acquires a large vev 
needn't be the object that gives the down-type fermions their masses. This allows a great 
deal of freedom in reconsidering the dynamics of the light fields in the theory (i.e., the 
charged Higgs) and whether they carry additional quantum numbers under some new group 
Gg. In particular, the charges of are only tied to the charges of S, not to the charges of 
down-type fermions. 

Secondly, the inclusion of the sister Higgs allows the presence of a new R-parity violat- 
ing operator, namely H^Tj^E, and with it comes a great deal of new phenomenology. In 
the presence of R-parity conservation, new LSPs are possible, including LSPs that carry a 
conserved charge and thus can potentially provide a dark matter (DM) candidate. 

Finally, the presence of these fields motivates the presence of additional down-type quarks, 
also charged under Gg- These G-quarks can have important consequences for signals such as 
gluino decays, and other opportunities for R-parity violation. Unlike usual 4th generations, 
these G-quarks decay via scalar (often Higgs) emission, leading to b-rich final states. 

Any one of these points motivates a consideration of these theories with extended Higgs 
sectors. We describe the construction of a sister MSSM in more detail below. In the rest of 
the paper we discuss the extension to a gauged sister model, before addressing each of the 
above mentioned points in more detail. 

A. A Sister Higgs 

A sister Higgs is a Higgs field that participates in EWSB but has no tree level couplings 
to matter. This is a natural outcome when the sister Higgs transforms as a non-trivial 
representation of some new symmetry Gg under which the SM fields are trivial. This may be 
a gauged or global symmetry. As we shall see shortly, although in its simplest manifestation 
Gs is not gauged, it is likely most interesting in a case when Gg is a gauge symmetry. 

To make this more concrete: we extend the MSSM by introducing two sister Higgs fields. 
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SU{3) X SU{2) X U{1) 



Y 



Gs [SUi2), 



S„ (1,2,1/2) r [2] 

Dg (3,1,-1/3) r[2] 

(1,2,-1/2) r[2] 

(3,1,1/3) r[2] 

$ (1,1,0) r[2] 

$ (1,1,0) r[2] 

Hu (1,2,1/2) 1 [1] 

Hd (1,2,-1/2) 1 [1] 

TABLE I: Charges of fields under the SM and the Sister group, Gs, for the generic case [the case 
of Gs = SU{2)]. Together and Dg (S^ and Dg) make a 5 (5) under SU{5). Note the usual 
MSSM fields Q, ^7^ L, and E"" are singlets under Gs- 

Su and Sd, which carry the same charge under the SM as Hu and Hd respectively, while also 
being vector-like under the new sister group, Gs- At this stage we do not specify whether 
Gs is a global or gauged symmetry. In order to allow couplings between the sister fields and 
the MSSM we also introduce fields that only carry charge under Gg- Furthermore, if these 
fields acquire a non-zero vacuum expectation value (vev) then there will be mass mixing 
between the sister Higgses and the Higgses of the MSSM. We denote these Gg breaking 
fields as $($), which are in the same representation of Gs as T,u{^d)- In addition to these 
fields, as we will see, it is natural to expect vector-like fields that transform as right-handed 
quarks and are also charged under Gg, which we dub G-quarks, Dg.Dg. Later, in Section |TT| 
we will investigate a particularly interesting case where the sister group is a gauged SU(2). 
The charges of all fields in both the general and this special case are shown in Table [Tj With 
this field content and assuming r 7^ f the simplest extension of the MSSM superpotential is 

+\u^Hj:d + Xd^^uHd 

+YuU'HuQ + YdD'HdQ + Y.E'HdL . (1) 



To make this completely general, we can add a term ^DgD'^ that we shall discuss shortly. 
Notice that if the representation r is different from the conjugate f, this superpotential 
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contains an additional U{1) symmetry, which needn't be accidental, which we denote f/(l)s. 
Under this symmetry S^, $ have charges 1/2 and S^, $ have charges —1/2. In the interesting 
case where Gs = SU{2), r is pseudo-real, and operators ^T^uHd + ^T^^Hu are also allowed, 
for instance, which break the ?7(l)s- Although the concept of a Sister Higgs is independent 
of Gs, if this group is too large (if the representations of S are larger than 4 dimensional) 
then the SM gauge couplings will reach a Landau pole below the GUT scale. 

The most natural origin for SM fermion masses is through the usual Yukawa couplings 
to Hu and as it is in the MSSM, and we have illustrated this in ([T]). However, with 
the inclusion of new fields acquiring vevs, it is possible that down-type fermion masses can 
instead come from the higher dimension operator ^HdQD^. This could be generated by 
integrating out a massive 5-1-5 pair, which are uncharged under G^, and have couplings 
EQV^ + 5$Srf. We will not focus on this case, but it may be natural to consider extending 
this model by the presence of additional 5 -|- 5's. 

II. GAUGED SISTERS AND THE HIGGS MASS 

As we have discussed, there has been great activity attempting to accommodate a Higgs 
at 125 GeV in supersymmetric theories. Such a large Higgs mass in the MSSM requires 
squark masses in the multi-TeV range and large A-terms, well beyond our expectations 
from naturalness [T2] . 

Various extensions of the MSSM attempt to alleviate this problem, as we have discussed. 
While the simplest possibility may simply be of a tuning of the MSSM, the simplest po- 
tentially natural model is likely the NMSSM. As is well known |8l [13], the RG flow of the 
SHuHd coupling in the NMSSM tends to make it run quite small in the IR. Unlike the top 
quark, which has SU (3) color contributing to its anomalous dimension, there is no strong 
interaction to drive this coupling large. Consequently, a 125 GeV Higgs lies near the edge 
of perturbativity without appealing to sizable radiative corrections from top loops. This 
leads us naturally to the possibility of gauging the sister group. If the sister group runs 
strong then the and couplings, which play the same role as SHuH^ in the NMSSM, 
can remain large and give a sizeable correction to the Higgs self coupling and therefore the 
Higgs mass. However, such an approach has a "why-now" problem: the coincidence of the 
strong coupling scale of the new gauge group and the weak scale. If instead the running 
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of the gauge coupling of the sister group, gs, were approximately conformal this why-now 
problem would be removed. An example of such a setup - where the beta function is zero 
at one loop - is = SU(2)s with^ six vector-like families of matter. 

With Gs = SU{2)s and $, $, S„, all in fundamental representations of SU{2)s a one- 
loop conformal beta function occurs if there are 3 more (□, □) pairs, the G-quarks. This 
is exactly the field content shown in Table [T] and is the same as what would be expected 
if the sister Higgses were embedded in a (5,5) of SU{5). Remarkably, simply by insisting 
on conformality as a solution to the "why-now" problem, we are naturally pushed into 
a unified (i.e. GUT) setup with the G-quarks. While the gauge coupling for SU{2)s is 
conformal at one-loop, the SM values are modified by the addition of two additional flavors, 
i.e., /33 = -l,/32 = 3,/3i=43/5. 

The low energy phenomenology of the sister Higgs scenario can be much more complicated 
than the usual MSSM (or NMSSM). In general, there are as many as six "Higgs" fields 
(i.e., SU{2)'[Y (g) U{1)y doublets) as well as four SM singlet scalars, in addition to a Z' 
and the associated fermionic fields. However, for the most part, we can understand the 
phenomenology by taking a simplified limit of the theory. As we will see this limit reduces 
the theory to a Type I 2HDM ^ 

We can take this limit as follows: let us begin by considering the superpotential in ([T]). 
This superpotential retains a global f/(l)s = f^(l)s + symmetry after SU{2)s breaking 
and an associated "sister charge". Within the SU{2)s doublet S^, for instance, only one 
component mixes with the Higgs, while the other component, carrying the sister charge, will 
not. When writing the scalar potential, we will designate with the * to indicate when 
the field carries sister charge, and similarly with the $, $ and S^j. In discussing the Higgs 
sector, we will for the moment assume the presence of this symmetry as it simplifies the 
physics by reducing by a factor of two the number of fields that mix with the Higgs boson. 

Since will take the role of the Ha for the purposes of raising the Higgs mass, is 
not important for the Higgs mass phenomenology, so for simplicity we set it to zero. We 

^ This group is also small enough that the additional matter charged under the SM does not lead to Landau 

poles in the SM gauge couplings. 
^ In a Type I 2HDM, there is a standard model like Higgs Hi which couples to fermions and a second Higgs 

H2 which does not, although it can acquire a vev. A mixing between them induces couplings for both 

mass eigenstates. 
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take SU{2)s to be broken by a = I'^cos^s and ($) = v<^sin^s. In general, we will be 
taking tan /^^ <^ 1, so the is acquiring the larger vev. For the moment, we then set — 0. 
We will assume that the cj) fields have large enough SUSY breaking masses that they can 
be integrated out without disturbing the phenomenology (i.e., the mass spectrum and the 
A^J/f-uErfp quartic) and set them to their vevs. E„ is for the most part a spectator, so we 
integrate it out with a SUSY breaking mass. 

The resulting theory, with these assumptions has a vev-acquiring sector composed of Ej^, 
Hd and H^, with the scalar potential, 

V = m^J/// + m^Ji/„r + (ml, + ^)|E/ + (ml, -^)|E^r (2) 
+ \Xuf\H^j:df + {BfiHuHd + Klil{^*HJla + h.c.) + ^i^^g^M + \dI ■ 

Here, we have taken the fermions to couple to the scalar sector through the usual Yukawa 

matrices YuQlf^Hu and YdQD'^Hd and Dmssm are the usual D-term contributions from 
SU{2)w and t/(l)y, and Dg are the residual SU{2)s D-term quartics after integrating out 
the 0, (f). 

We can further integrate out in the large m\d limit by — BjiH'^/m^jd^ which yields 
the final simplified setup. 

2 2 

V = m|Ji/„r + (^L + ^)l^<^r + (^L - "^Wd? (3) 
+\K\^\HJld? + (A„/x;($)*i/„Erf + h.c.) + \dI,ssm + \dI , 

where the fermion masses now arise in the effective theory via the terms YuHuQV^ and 
Yd^^H*QD'^. Put simply, the low energy effective theory is a Type-I 2HDM with an 
additional spectator E^. Although E^ does not mix with these fields nor acquires a vev, we 
have left it in because of the relationship between its mass and that of the U{l)s neutral 
E(i (which can lead to tensions in some regions of parameter space). Prom here, the scalar 
phenomenology is fairly simple. acquires a vev (Hy) — v sin /?, while (E,^) = v cos (3, and 
the Higgs quartic has an NMSSM hke contribution proportional to A^. 

The NMSSM-like boost to the Higgs mass will depend on how large A„ can be. The (3 
function for A„ at one loop is 

Pxu = A„ (^32/2 + 5X1 - lal - ^gl - ^gi^ , (4) 
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FIG. 1: IR values of tan (3 and Xu when scanning over the range [0,4] for the UV inputs gs, Xu, and 
yt- Different ranges for gs{lTeV) are denoted by different symbols, o: gg < 0.5, : 0.5 < gs < ^, 
1 ^ 55 < 1-5, and O: gs > 1.5. The lines show the parameters for which a higgs of 125 GeV 
is achieved at tree level (solid), with stops of 250 GeV (dashed), 500 GeV (dotted) and 1 TeV 
(dot-dashed) . 

where (71,2,3 are the standard model couplings, and Qs is the coupling of the new sister 
group. These expressions make clear the simple point: the presence of Qs acts to counter the 
tendency of A„ (and similarly for A^) to run small in the IR. A simple question, therefore, is 
how large can the Higgs mass be in this theory, consistent with perturbativity up to the GUT 
scale? Because some couplings are sizeable, it is best to calculate the RG flow to low energies 
at two-loop. To do this, we use the publicly available SARAH package [I5l [16], to evolve the 
couplings within our theory from Mqut = 2 x 10^^ GeV to 1 TeV. We generate a large set of 
random high scale values of Qs, yt-, Xu, all lie in the range [0, 4], and run down to low energies. 
We define sin/3 = ^P^^^y where v = 174 GeV and we take mt(lTeV) = 150.7 GeV pjj. 
Points with sin /3 > 1 are disregarded as in conflict with observation. 

The results of our scan are shown in Figure [T} While one must be careful with this plot 
(as there is no sense of a measure on this scatter plot) we can get a sense of correlations 
between parameters. First, we see that it is quite easy to have very large values of A„ 
at low energies. Taking ag^McuT) ^ 0.18, i.e., gsiMcur) > 1-5, we can have A^ ~ 1.25 
at low energies. This is sufficiently large that no additional correction from stop loops is 
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FIG. 2: Example spectrum of fermions, in our simplified limit. States in red are charged under the 
global U{l)s- No special ordering of mass scales is implied. 

necessary. However, including even moderate stop contributions is easy. With as = .02 (i.e., 
comparable to the SU{2)w coupling) A„ is easily large enough to get mn = 125 GeV with 
only ~ 500 GeV stops. 

A final point on which we make a small comment, is that $ plays the role of the usual 5" 
in the NMSSM and, as recently pointed out by [3J, for large Xu the theory can become more 
natural, in that the soft Higgs mass can be larger, and then mixing the SM singlet scalar 
(in this case (p) with the Higgs, its physical mass can be lowered at the expense of just a 
moderate tuning. Although we do not explore this quantitatively here, it is clear that the 
ingredients (a large A„ and a neutral scalar mixing with the Higgs) are present. 



A. New Phenomenology in the Simplified Limit 

Even in this limit where we have decoupled 0,0 and there are very important differ- 
ences from the MSSM. First, because the second Higgs field in the low energy theory, S^, 
has no direct couplings to SM fermions, the mass constraints that apply to in Type-H 
2HDM do not directly apply. In particular, the charged Higgs can be < 150GeV without 
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FIG. 3: Example spectrum of scalars and vector bosons, in our simplified limit. States in red are 
charged under the global C/(l)s. 



necessarily contributing too largely to the decay of the top, nor excessively to 6 — )■ 57 |18j . 

However, unlike a pure Type-I 2HDM, with Gs = SU{2), the sister Higgs is tied to a 
second field, namely its sister-charged partner E^. In particular, the mass of the charged 
Higgs if^ (which is dominantly S^) is related to the mass of the sister-charged by 



m^± — m'^ ^± ~ rri^z (assuming ($) ^ ($))• With larger ($), this becomes Tn'jj± — 



2 

m ± 



(cos2/3s + 2A^ sin^ /Ss/fiff ) "^1^ (assuming tan/3 ^ 1). Either way, for heavy enough Zs and 
light enough if^, this will induce a charge-breaking negative squared mass to T.'^^ . 

Probably the simplest way out is to break the global sister charge f/(l)s by including an 
operator such as XZ^HuTia- The F-term for Hu then gives an additional contribution to the 
mass of S^, avoiding the unwanted charge-breaking negative ^± . This operator, combined 
with A„$if„Erf, will also contribute to the squared-mass of through the F-term for S^. 
In order to avoid fine-tuning, ($) and ($) should not be too large, implying that Zg should 
not be too heavy relative to the EW scale. Moreover, including this sister-charge breaking 
term will affect the running of A^, but given that large values of A^ are already generally 
quite natural, this is unlikely to be problematic. 

Another key element of this setup is the presence of new gauge forces, Z'^ and W^. Because 
Erf is charged under both SU (2)^ and 5'[/(2)w, its vev will induce mixing between Z° and the 
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FIG. 4: Bounds on Zg — Z mixing expressed in terms of mz^, tan/3 and Qs, from precision elec- 
troweak constraints of the p-parameter (red) and direct collider searches for dilepton resonances 
|19j (blue). The later are derived assuming that Zg has the same branching ratios to leptons as 
the SM Z^. The region below each line is excluded. 

SM Z". Such mixing is constrained by precision electroweak (PEW) measurements of the p- 
parameter (although, as has been pointed out previously, e.g., by [20J, this can actually help 
agreement with the standard model) as well as collider searches for Z' resonances. ATLAS 
performed a search for high-mass dilepton resonances in [19] . We can conservatively re-cast 
the results of that search in terms of Mz^, Qs and tan/3, assuming that Zg decays as a 
sequential Z' {i.e. assuming that its branching ratios to SM fermions is identical to that 
of the SM Z°). These bounds, as well as the PEW constraints, are shown in Figjij Note 
that even for large Qg and moderate tan /3, Zg can be relatively light. Finally, we point out 
that the direct bounds from dilepton resonances are significantly weakened if the dominant 
decay mode of Zg is Zg — )■ Z^h^,W^H^ or even (f)*(j), which is naturally the case if such 
on-shell final states are kinematically allowed and tan /3 is not too large. Despite this, the 
PEW constraints are stronger and hold regardless of how Zg decays. 

An interesting consequence of Zg — > being the dominant decay mode of Zg is that it 
provides a sizable contribution to Z^h^ production at the Tevatron. Saturating the precision 
electroweak constraint on Z — Zg mixing, the Zg production can be comparable to that of 
SM Z^h^ production at the Tevatron for relatively light Zg, as shown in Figjsj At the LHC 
such effect is much smaller, however. 
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FIG. 5: Saturating the bounds for the p-parameter, there can be a sizeable contribution to Z^h^ 
production at the Tevatron through resonant Zg production. At the LHC, it is a small correction 
to the SM process. 

Finally, we comment on the implications of the Sister Higgs framework on the phe- 
nomenology of dark matter. As was pointed out in [21], in the phenomeno logically viable 
region of MSSM parameter space, achieving the correct relic abundance for LSP dark mat- 
ter requires considerably tuning between parameters. There are several viable regions where 
the so-called "well tempered" neutralino lives, with just the right admixture of bino/wino 
or bino/higgsino to allow the LSP to be light (m^ ~ (9(100) GeV) and have the correct relic 
abundance. In the case of the Sister MSSM there are many new neutralino states that can 
mix with those of the MSSM; in general the neutralino mass matrix is now 13 x 13. Many 
of these new states carry charge under Gg and have new interactions, either through the 
superpotential or, if Gg is gauged, Zg/Wg exchange. The new gauge bosons can be light 
and have sizeable gauge couplings (Fig. and so have considerable effect on the process 
of thermal DM freezeout. The possibility of Dirac dark matter charged under the residual 
sister U{l)g will be explored at length shortly, our point here is simply that the neutralino 
sector is much richer, which allows for a broad range of possibilities for DM. 

III. R-PARITY VIOLATION AND A SISTER HIGGS 

With the addition of the new fields charged under Gg, the symmetry structure quickly 
grows complex. As a consequence, some of the new fields may be stable. At the same time, 
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TABLE II: Charges of fields under the global symmetries. The non-MSSM fields are also charged 
under a gauged Gg- 

we shall see that an intriguing new opportunity for R-parity violation appears in this theory. 
Let us begin by considering only the superpotential of eq.([T]). This superpotential has the 
gauge symmetry group SM ® SU (2)^, but has a global symmetry U{1)b® U{1)l ® U{1)y. ® 
U{1)g X Pr- B and L are the usual baryon and lepton numbers, under which only the matter 
fields of the MSSM are charged, the new fields charged under SU (2)s are taken to be neutral. 
G is a conserved quantum number for the G-quarks, Pr is the usual R-parity. With this 
superpotential, $, and can be all even or all odd under Pr, and similarly the parity 
of Dg is unassigned. However, to avoid spontaneous breaking of R-parity, the components 
of if, S, and $ scalars which acquire a vev must be even under R-parity. Finally, is 
a global symmetry under which $, have charges —1/2, while $, have charges +1/2. 

When SU{2)s is broken by $ and/or $ vevs, the SU{2)s (g) f^(l)s is broken to a global 
U{l)s, which we refer to as a sister charge. Under this residual symmetry, one component of 
$ is charged, while the component getting a vev is neutral, and similarly for $. Likewise, one 
SU{2)l doublet of mixes with Hfi and is neutral under U{l)s, while another is charged 
and does not mix, and again with and H^. Thus, in this theory, with no additional 
operators, the hghtest baryon, lepton, R-parity odd field, sister-charged field and G-quark 
are all stable. The charges of the fields are summarized in Table |TTj 

The presence of additional operators can break the symmetries further. In particular it 
is possible to break R-parity in several ways. Below we discuss a particularly interesting 
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case, with novel phenomenology, that breaks R-parity in a new way in the lepton sector, 
but before doing so we discuss an alternative approach, with less novel phenomenology, of 
R-parity breaking in the hadron sector. 



A. Hadronic RPV 

It is simple to break U{1)g, for instance by including a small mass-mixing term^ 

^D^D". (5) 

This operator has charge under B, G and S, and its inclusion breaks one linear combination 
of these symmetries. The two unbroken symmetries can be taken to be the combinations 
B' = B + G/3 and E' = E + G'/2. ^ Under this modified baryon symmetry, B', the field Dg 
carries charge 1/3. This will allow G-quarks to decay, assuming that there is some lighter 
particle with the residual U{l)s charge (such as S or $). 

Another operator which affects the U{1)g is DgD^W^. This operator is charged under 
U{1)g and U{1)b, and thus breaks a combination. We can think of the remaining symmetry 
as a modified baryon number, B" = B — G/6, with this assignment Dg has new baryon 
number —1/6. Note that the G-quarks cannot decay with just this operator, since there 
are no states lighter than them with B" = ±1/6. Equivalently, it retains a G-quark parity. 
Including this operator together with ([s]) breaks U{1)g ® f^(l)s ® ^^(1)^ down to a single 
f/(l). We can think of the residual number as a baryon number, under which D'^, Dg, $ have 
charges of —1/3, —1/6 and —1/2. 

When $ acquires a vev it breaks SU{2)s f^(l)s down to U{l)s, which we can now 
identify as baryon number, B'" = B" — (T^ + S'). That is, different components of Dg under 
the SU{2)s transformation have different baryon numbers. Simply put, since the operator 
is Dg^Dg^W^ (with 1,2 indexing the SU{2)s component), if we set D^^ as the field that mass 
mixes, then we can treat D'f as a diquark, and the residual sister symmetry as baryon 



number. Since, in a similar fashion to the lepton RPV of Section III we can redefine baryon 



^ We include this rather than ^DgD'^ at least initially, as Dg is considered to be in a multiplet with and 
thus to have the same Gs representation. Assuming the presence of even an approximate ^DgD'^ 
should be more significant than ^DgD'^. 

^ Note that, alternatively, we could have included the operator QD'^^^, which would have the same effect, 
although it breaks the three charges down to a different linear combination. 
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number and therefore R-parity to include the sister charge, this operator does not allow the 
LSP to decay as a consequence. 

We can achieve R-parity violation in this fashion if we include the operator ^DgD'^, along 
with '^DgD'^ + DgBgW^, which explicitly breaks U{l)j]. In essence, this mass mixes both 
components of the G-quarks, so that DgBglf^, in the mass eigenstate basis, yields the usual 
UDD RPV operator. Alternatively, if there is an additional 5 + 5 uncharged under SU{2)s, 
we could integrate it out and leave the operator 

^ ^DID^U^ (6) 



Mr-' 3 



'55 

In the presence of (|5]), this acts as the usual baryon number violating operator UDD. 
Interestingly, this preserves f/(l)s and thus the lightest sister particle (or LSiP) remains a 
potential dark matter candidate. 



Finally, we comment that including both terms such as ErfS^i^"^, discussed in section III B 



as well as ^DgD'^ + D'gD'gU^ can lead to B + L violating (but B — L conserving and i?-parity 
conserving) operators such as U D'^D'gEHd^d- Since these operators preserve B — L, they 
will not lead to proton decay, but could lead to neutron decay inside a nucleus. While 
the constraints on such operators are weaker than with B — L violation, they are clearly 
constrained, and we leaved detailed studies for this hybrid scenario for future study. 



B. Leptonic RPV 

Within the MSSM, the operator H^HdE'^ is never discussed, and for good reason - since 
the SU{2)w indices of Hd are contracted antisymmetrically, it vanishes identically. However, 
as soon as there is an additional flavor of Hd, the operator H\H^E'^ can be present. This 
sister-Higgs R-parity violation is qualitatively different from usual RPV in the MSSM, as 
we shall discuss. Let us first see how this arises in the more complicated scenario where the 
sister Higgs is an SU{2)s doublet. 

We can naturally include the superpotential operator 

i^i^d^dE'i . (7) 

This breaks ?7(l)s ® down to ?7(1)l' where L' = L — H and under which Hd carries 

lepton number 1/2. In addition this operator determines a unique R-parity assignment 
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within the S fields that leaves R-parity preserved. When $ acquires a vev, we can think of 
the residual sister charge as identified with lepton number. Hence, this operator does not 
break R-parity without some additional effect. More explicitly, there is a new, preserved, 
R-parity we can define, B! = (_i)3(-B-Qs)+2s^ where Qs = + L' and s is the spin of the 
particle. The simplest modification that leads to R-parity violation is to include the terms 
^HyT^d and/or in the superpotential. These explicitly break f/(l)s, which means 

there is no residual f/(l)s, and thus both components of will ultimately mix with Higgs 
bosons. There are no stable states in the theory that may act as dark matter. 

Alternatively, if again we consider a 5 + 5 uncharged under SU{2)s, we could integrate 
it out and leave the operator 

-^^^^HaEt . (8) 

This operator preserves f/(l)s, and hence ?7(l)s, but it is not possible to define a preserved 
R-parity in the presence of ([T]) and ([s]). This operator breaks R-parity (and lepton number), 
and thus we are left with an interesting possibility: that R-parity is broken, but that a 
residual global symmetry is intact, leaving a potential candidate for dark matter. In this 
case the collider production of DM will either be through electroweak processes, or through 
the decay of the G-quarks which have QCD production, but may be very heavy and so 
kinematically suppressed. 

In either of these two cases, stable DM candidate or not, there are interesting implications 
for collider phenomenology. R-parity violation (RPV) will remove the usual ^-t signature 
of SUSY, but, unlike hadronic RPV, where the $x is converted to jets (and can thus be 
extremely challenging to discover), here, the $j' is converted to leptons+jets, which is much 
more tractable. Furthermore, under the assumption that some fiavor structure dictates the 
sizes of couplings, the final state lepton is essentially all r. Then the final state will contain 
non-negligible and the SUSY signatures, while more challenging, are not necessarily 
impossible. If the dominant decays are to e, /i then the requirement of many new 
physics searches may not be met and the constraints will be weaker. We discuss below the 
possible signatures for decays for various LSPs in the presence of this form of leptonic RPV, 
and outline which existing searches may be sensitive to these signatures. We leave a detailed 
study of efficiencies and resulting bounds for future work [22j . 



16 




FIG. 6: Squark and stop decay modes in the leptonic sister RPV scenario. 
C. Constraints on Sister RPV 

Both ([7]) and ^ will lead to mixing between the RH leptons and the Higgsinos. We 
concentrate for concreteness on the effects of ([t]) with one or both of + ^Hylld- 

One linear combination of leptons, the combination ^ will have a mass contribution 

nfv"^ cos^ /3/My,. Unlike the conventional LH^ and LLE'^ RPV operators these contribu- 
tions do not affect the left-handed leptons, meaning that constraints from neutrino mass 
are much weaker. Likewise, since there is only a single lepton appearing in this operator, 
constraints on dilepton resonances are not expected to be relevant. ^ 

Rare decays are expected, but because this operator couples to right-handed leptons only, 
the transition will have a chiral suppression, i.e., we expect a dipole operator for — )■ £^7 
{li = fi,T, if = /i,e) of Dif ~ miKiKfe/16n^M^. Current bounds on the lepton flavor 
violating branching ratios place constraints on the Ki of, 

^ and < 2 X 10"^ ^ . (9) 

looGevy \ T i\ ^ ViooGevy ^ ' 

Squark LSP 

A squark LSP in the presence of sister RPV is an interesting case to consider. Given 
that the operator H\H^Ef mixes the chargino with a charged lepton, the dominant decay 
mode of a squark LSP may be g — )■ q'i^ ( Fig. |6|, similar to the decay induced by the 
RPV operator QLD^ in the MSSM. Note, however, that contrary to QLD'^, the sister RPV 
operator H^HjEf does not induce the decay q qu^ at tree-level. 

If the squarks belong to the first or second generation and dominantly decay to g — ?■ q'e^ 
or g — )■ g'/i^5 straightforward bounds on rriq can be directly extracted from leptoquark 

^ The exception is if the LSP is a sneutrino, which can decay to two leptons. 
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searches. ATLAS carried out searches for first and second generation leptoquarks with 1 fb~ 
of luminosity [231 El], and constrained the masses of leptoquarks to he itilq > 685 GeV for 
LQ fij and rriLQ > 660 GeV for LQ — )■ ej. CMS also searched for second generation 
leptoquarks with 2 fb^^ of luminosity [25], and placed the bound niLQ > 632 GeV for 
LQ fij. 

The bounds above on leptoquark masses directly apply to squarks decaying to those same 
final states, since the production cross section for squarks and leptoquarks are effectively the 
same in the limit in which the gluino is heavy and has a negligible contribution to squark pair 
production. However, in generic scenarios in which the first and second generation of right 
and left handed squarks are degenerate, the squark pair production cross section is a factor of 
8 larger than leptoquark pair production. It is straightforward to extract the limit on squark 
masses in this case: the ATLAS search bounds the production cross section of g — )■ q'e^ to 
be a{pp — > qq* — > e"'"e~jj) < 5 x 10^^ pb, implying that rriq > 880 GeV. Similarly, the 
strongest constraint on g — )■ q'fi^ comes from CMS. It bounds the production cross section 
to be a{pp — > qq* — > fi'^fJ^'jj) < 3 x 10^^ pb, which translates into rriq > 900 GeV. 

Those searches cannot, however, be as straightforwardly applied to the case in which the 
dominant squark decay mode is g — )■ q'r^. In that case searches specifically targeting new 
physics signatures with jets and two tau leptons place stronger bounds. ATLAS performed 
a search for events with jets, at least two tau leptons and missing energy using 2 fb~^ of data 
[26] . We re-interpreted the results of those searches in order to place bounds on first and 
second generation squarks decaying as g — )► q'r"^, as well as stop LSPs decaying as t — )■ 6r+. 
We used MadGraphS to simulate pp — )■ gg* and pp — )■ it* at the 7 TeV LHC. Those events 
were passed through Pythia6 for decay, showering and hadronization, and then handed to 
PGS4 for a rough simulation of the detector response. 

The ATLAS search in [26] selected events with at least two jets, the leading jet with 
transverse momentum pt > 130 GeV and the sub-leading jet with transverse momentum 
Pt > 30 GeV. In addition, it required two hadronic tau candidates with pt > 20 GeV and 
missing transverse energy > 130 GeV. It also placed a cut A(f){$rp,pTi^2) > 0.4 between 
the missing transverse momentum and the two leading jets, and rrics > 700 GeV. In Figure [7] 
we display the distribution of events in terms of the sum of the transverse mass of the two tau 
candidates, for a few benchmark signals. The signal region chosen by ATLAS corresponds 
to + m'!^ > 80 GeV, and it bounds the number of signal events in that region to be 
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FIG. 7: Distribution of + m!^ for the event selection of the ATLAS search in [26]. The SM 
background, as estimated in [26j, is displayed in gray, and the measured data is shown as black 
dots. Two benchmark signals with masses at the boundary of the this analysis' reach are displayed 
as well: (left) degenerate 1st and 2nd generation squark LSP's decaying as g — t- q'r^, and (right) 
stop LSP decaying as t — 6r+. 

< 6 events. For the scenario of degenerate first and second generation squarks decaying as 
q — )■ q'r"^, that excludes squark masses nig < 720 GeV. For the stop scenario with decay 
t — ?► 6r+ the exclusion is < 540 GeV. This bound is surprisingly similar to the one 
obtained by the CMS dedicated search to this signature [23, which constrains stop masses 
to be mi < 525 GeV. Searches for charged higges from stop decays, t — )■ bH~^ — )■ 6(r"'"i/), 
such as [28] , also yield non-trivial bounds. We estimate that the bound from this particular 
search to be < 300 GeV, hence weaker than the bound from [26l |27] . 



Neutralino LSP 

In the presence of the leptonic sister RPV operator H^H^Ef, a neutralino LSP can decay 
to i^H^. As we discussed before, the phenomenology of the light Higgs sector is that of 
a type-I 2IIDM, and therefore will decay dominantly to tb (or Wbb), tv and cs, with 
branching ratios determined by its mass and whether it is allowed to be on-shell or not. If 
is too heavy, the decay Xi ~^ d-^W^ might dominate. Either way, the expected signatures 
from neutralinos at the bottom of SUSY cascades should contain a large number of leptons. 
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b,s,T q,t 
FIG. 8: Neutralino LSP decay modes induced by leptonic sister RPV. 




FIG. 9: Chargino LSP decay modes induced by leptonic sister RPV. 

If these neutralinos are being produced through strong interactions (such as gluinos or 
squarks), these scenarios should be significantly constrained by same-sign dileptons [291433] 
and multilepton searches [31H39] . which have very low backgrounds and hence have very 
powerful sensitivity even for models with very low branching ratio to multileptons. Other 
potentially constraining searches are same-sign tops (if the decay Xi ^ t^H^ — >■ i^(tb)^ is 
dominant), first and second generation leptoquarks [23ti25] . or GMSB with LEV |40j . 

Chargino LSP 

It is also interesting to consider the situation where the chargino is the LSP (see [H] for a 
discussion of scenarios), or when its decay to a neutralino is suppressed due to phase space. 
In this case, the decay — ^ ^^^S ^ill be quite challenging, see Figure |9} If hj. couples 
to up quarks, it is possible that we would get a sizeable signal of WW (through its mixing 
with the Higgs), but if unavailable, then cc would dominate. Or, if coupling through down 
quarks, bb, which can perhaps aid by adding b-tags along with hard leptons (possibly rs). 
Alternative decays to W^u or i^Z^ should be constrained constrained from standard SUSY 
searches in jets and/or multileptons. 
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Slepton/Sneutrino LSP 



Slepton or sneutrino LSPs decay in a similar fashion as in standard leptonic RPV induced 

uction have been recently 
^z/ (where ^^ = or /x^) 



by the LLE operator, see Fig, 10 Searches for direct slepton production have been recently 
released by ATLAS pi2| and constrain sleptons decaying as £^ — )■ 
to be mi± > 185 GeV. These searches, however, cannot constrain sneutrino pair-production 



due to the absence of missing energy. However, multi-lepton searches requiring 4 leptons 
should have sensitivity to — )■ 4i. 





FIG. 10: Slepton and Sneutrino LSP decays induced by leptonic sister RPV 



Gluino LSP 

Gluino LSPs will decay to dijets plus a charged lepton. If the charged lepton is a r^, 
standard jets+-^T plus 0, 1 or 2 charged leptons should have sensitivity to pair-produced 
gluinos. For gluinos decaying to jets plus or fi^, however, the absence of missing energy 
makes this signature more challenging. Limits from current searches on leptoquarks, black- 
holes or same-sign dileptons may be re-cast for gg — > {qq(-^){qq('-^)] however dedicated 
searches may have a potentially much higher reach. 

IV. DARK MATTER 

One appealing element of SUSY is the presence of a dark matter candidate in the form 
of the LSP. Under the assumption of R-parity conservation, the neutralino often serves as 
a WIMP candidate. In this theory, because of the sister Higgsinos, the state forming the 
WIMP is more complicated, with the new forces altering the DM freezeout story. While a 
thorough study of dark matter is beyond our scope, some discussion is clearly worthwhile. 
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A. Dark Matter with R-Parity Violation 



In general, there is a tradeoff in models between the presence of RPV (which can be an 
explanation of the so-far elusive nature of supersymmetry) and a dark matter candidate. 
The standard lore is that since the LSP is the dark matter candidate, violating R-parity 
(and thus removing the $rp signal) will remove a cosmologically stable particle from the 
theory. However, within the sister MSSM, that lore is easily violated by the presence of a 
residual global symmetry. 

In the Sister Higgs scenario, that symmetry can be the global f/(l)s, which, if preserved, 
makes the lightest state charged under this symmetry (the "LSiP") stable and hence a 
potential dark matter candidate. If the U{l)s is broken down to a Z2, then the LSiP is 
Majorana. As we shall see, the Dirac case is quite constrained by direct detection, but the 
Majorana case is much less so. 

Among the fermions carrying f/(l)s charge, there are three Dirac states comprised of 
0g = {^1,^2)1 = (E„i,Ed2), and Wg = (Ws+,Ws-). The masses and mixings of those 
fields can be extracted from: 

C D [$1 S„i Ws+] -Kvsp ^,($2) 

/is Srf2 • (10) 

gsi^i) QsVCp -M, 

Those fields couple vectorially to the Standard Model boson, either directly though 
their S^- components, or through — Z^ mixing. If the lightest mass eigenstate is the 
dark matter, its coupling to the Z^ is constrained by direct detection experiments. The 
LSiP-nwc/eon cross section induced by Z^ t-channel exchange is: 
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{A-Zf 



(11) 



where klsip-z is the dark matter coupling to /i„ is the dark matter-nucleon reduced 
mass and A (Z) is the nucleus atomic mass (number). This cross section is constrained by 
XENONIOO pSj to be cr„ < 10~^^cm2. This translates into a bound on the dark matter- Z° 
coupling of KLSiP-z ^ 3.7 x 10""^, which rules out as a Dirac dark matter candidate since 
it is charged under SU(2)w and its coupling to Z^ {f^ts-z" — 92/'^Cw) violates the direct 
detection bound by three orders of magnitude. This bound also severely constrains 0^ as a 
Dirac dark matter candidate, since it mixes significantly with S., through Au(if„)Srf$. 
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These bounds can be evaded however, by adding small Majorana mass terms, for instance: 

£ D 5mi<li<li + (5m2$25>2, (12) 

or analogous terms for S^. This term breaks U{l)s but preserves a Z2 parity that stabilizes 
the dark matter. It also splits the mass eigenstates by 6mi + 6m2- Since those states are 
Majorana, their Z-mediated elastic scattering is spin-dependent, and moreover suppressed 
by 6m/ fi^. The spin- in dependent scattering is inelastic and hence can be suppressed or 
even shut off for large enough mass splittings Note that it is non-trivial to achieve 

this term as superpotential couplings -|- <l'<l> vanish identically by the antisymmetry of 
the SU{2)s indices. Thus, generating this term would naturally involve SUSY breaking or 
higher dimension operators. We shall not pursue this except to note that this is not a simple 
model addition. 

In the absence of such [/(l)s-breaking majorana masses, however, the only viable Dirac 
dark matter is a state that is almost purely sister gaugino, Wg. The direct detection bound, 
in this case, translates into constraints on: 

• the mixing of Wg with E^: 

^_9l^_k^:,±l3^^S.7xlO-\ (13) 

2 Cyj 2 

• the mixing of the Standard Model with the sister gauge boson Z^: 

gsxOz-z, = / . . 2 ^ 2 cQs'/^ < 3.7 xlO-^ (14) 



The later constraint is illustrated in the left-hand plot of Fig, 11 for tan /3 = 4. In the gray 



colored region the Z — Zg mixing violates the bound in (14), predicting a direct detection 



rate larger than presently allowed. Constraint (13) on the mixing between Wg and Hg can 
be satisfied with large /i^ and /is terms (see eq.(lO)). Fig, 11 illustrates the allowed regions 
in the mz^ — Qg parameter space for /i<^ = /is = 0.9 TeV, 1 TeV and 1.2 TeV. One can see 
that the larger and /is are, the more 0^, are decoupled, suppressing their mixing with 
Wg and allowing a larger parameter space in mz^ — Qg- Thus, insisting on a dark matter 
candidate of this type pushes us to a very unnatural corner of parameter space, since /i<^ and 
/is are directly tied to electroweak symmetry breaking, and their TeV scale values imply 
large cancellations in order to get the correct Z^ mass. 
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FIG. 11: (Left) Direct detection bounds on gg vs mz^ fixing tan/3=4 and Au=0.6, assuming that 
dark matter is dominantly tlie gaugino LSiP, Ws- The gray region is excluded due to large Z — Zg 
mixing, and the colored regions are allowed when the mixing of Ws with is sufficiently suppressed, 
which can be achieved with ~ C'(l) TeV. (Right) Wg relic abundance, assuming mz^ = 

600 GeV. 

In order to complete our dark matter discussion, we also consider the dark matter relic 
abundance. The dominant WgWs annihilation channel will depend on whether there are 
states with SU{2)s interactions that are lighter than Ws, such as the scalar components of $ 
or $, or G-quarks. If annihilation into those (on-shell) states is kinematically allowed, they 
will dominate due to large (7^- couplings. If not, then the dominant annihilation channel will 
be VFsPVs — i- — i- Z^h^. The final state W^H^ is kinematically forbidden since the mass 
of the charged Higgs, controlled by Xufi,j,{^), is at the TeV scale. Annihilation through an 
s-channel Z^ to Standard Model fermions is subdominant relative to the Z^h^ final state 
when Z — Zg mixing is suppressed. However, for mz^ ^ 350 GeV, the annihilation cross 
section to SM fermions is larger than 10% of the cross section to Z^h^ , so we include it as 
well in what follows. 
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FIG. 12: Relative IFsll^s annihilation cross section to Z^h^ vs. light SM fermions as a function of 
mzs- The LSiP mass is set to = mzj2. 



The annihilation cross section for ll^oW^, — > — i- Z'^h'^ is: 
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(15) 
(16) 



Ws Ws/ ^ Ws 

where as = gg/{4:n), a is the mixing angle that determines the mass eigenstates of the CP- 
even higges and H^, and we can ignore the Z^-width since, in the region of parameter 
space we are considering, it is < 0{MeV). Moreover, in the limit we are considering the 
CP-odd higgs is very heavy, so that sin a ^ — cos /3. Finally, 



^{x, y) = [1 + 2{5x -y) + {x- yf) ^/l + x"^ + y"^ - 2{x + y + xy). 
Annihilation into SM fermions through s-channel Z^ or Z^ is given by 



(17) 
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From (15) and (18), we see that the dominant annihilation channel {Z^h^ vs //) is 
determined primarily by the mass of Zg (there is a smaller dependence on m^y^ that is only 
important near threshold). Fig jT2] displays the relative fraction between the two annihilation 
channels (assuming m^y^ ~ mzs/2). As mentioned previously, for mzs ^ 350GeV, the 
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branching to // drops below 10%. However, it is the dominant annihilation channel for 
mz^ < 240GeV. 

In order to get the observed relic abundance the total annihilation cross section must be 



C(l) pb i.e. {av) = {crv)zOhO + W'^)ff ~ 1 pb- The right-hand plot of Fig, 11 displays the 
values of Qs vs. m^y^ for which the correct relic abundance is obtained, fixing mz^ = 600 GeV 
and tan /3 = 4. 

Alternatively, if we insist on naturalness and still wish to preserve the global U{l)s 
symmetry, we can give up on having the LSiP as dark matter and instead interpret the 
direct detection bounds as a constraint on the LSiP relic abundance. In that case the bound 
can be re-cast as: 

«LSiP-z- ^ (3.7 X 10-^)^ (19) 
Po 

where p is the local dark matter density and po ~ 0.3 GeV/cm^. For an LSiP with elec- 
troweak charge such as S.,, that implies: 

P < 10-' Po. (20) 

Achieving such an under-abundant relic density is challenging. As mentioned before, the 
annihilation cross section can be significantly increased by allowing the LSiP to go to on- 
shell states that are charged under SU{2)s, such as 0's or G-quarks. Even so, the increase 
is cross section relative to the Z'^h^ channel scales as cos~^ /3, a gain of ~ 3 x 10^ over what 



is shown in FigjTTj Although an additional boost may be achieved with large values of gs, 
naively resonant annihilation through would still be necessary in order to deplete the 
LSiP by six orders of magnitude relative to ^dm- A tension then appears since the Zg width 
is also increased to Tz, ~ 0{1 — 10) GeV, intrinsically limiting the efficiency of resonant 
annihilation. Other channels such as t-channel annihilation into ZgZg cannot easily achieve 
the required under- adundance either, even for as ~ 1. Note that need for QlsiP ^ 10~^ ^dm 
only holds if the LSiP is dominantly the SU{2)-ew doublet E^. Hence, stable Dirac LSiPs 
that are dominantly $s ot Wg may still be viable. 

In summary, there seems to be a basic tension to achieving a combination of a) having 
an intact U{l)s, b) having a stable Dirac particle charged under U{l)s and SU{2)ew, and 
c) naturalness. The easiest ways out are to break the f/(l)s completely or to a Z2 (in both 
cases the LSP will be Majorana) and making the spin-independent cross section inelastic, or 
to preserve the U{l)s but identify it with an existing SM charge such as U{1)b or U{1)l, in 
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which case any heavy particles could decay to SM quarks or leptons, so that direct detection 
constraints will not apply. 



B. Bosonic LSiP Dark Matter 



So far we have focused on fermionic LSiP dark matter, but it is worth commenting on 
bosonic LSiP dark matter as well. Many of the same issues will arise here, again due to 
Z — Z' mixing. 

The candidates for a bosonic LSiP WIMP are the '■^^ charged fields, i.e., 0*^**^ a^'^^ 

{s) is) 

and Gu among the scalars and Wp, among the vectors. Examining the spectrum of Figure 
|3]we see an immediate tension in our simplified limit. Although the usual SU{2) D-term 
pushes (Tq''^ down in mass, the NMSSM-like term raises the electrically neutral component 
in general above the electrically charged component. Thus, even to find a reasonable DM 
candidate we must explore a broader range of parameter space. 

Even so, for a suitable candidate, because it is charged under SU{2)s, the Z — Z' mixing 
will induce the same large Z-exchange direct detection process, meaning that only a very 
light WIMP would be allowed, or under a situation where f/(l)s is broken down to Z\. 

An interesting possibility in the broader space would be if the W^'^^ was the LSiP, yielding 
an interesting possibility of vector DM. It, too, would interact through the Z, and so similar 
model building constraints apply. 

Clearly the model space is broad, and is beyond the scope of what we can accomplish 
here. We have only begun the discussion of DM in Sister Higgs models. 



V. COLLIDER IMPLICATIONS OF AN LSiP LSP 



If the LSP carries sister charge, the collider phenomenology of pair produced superpart- 



ners is modified in an interesting way. Fig, 13 illustrates that for a squark decay. Since the 
squark does not carry sister-charge, it can only decay to final states with pairs of sister- 
charged particles. Spin-statistics requires that the final state contains the lightest sister 
fermion and the lightest sister boson (note that both are stable). If such decay is mediated 
by an off-shell neutralino, its kinematics will be that of a 3-body decay, where two of the 
three final states are "invisible". 
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An analogous situation occurs in the MSSM with an off-shell Bino mediating a 3-body 
decay of the squark to a quark, a neutrino and a sneutrino. Generically a sneutrino LSP 
would be severely ruled out by direct detection; however those bounds could be evaded in 
scenarios with mixed-sneutrino LSP [15], or a sneutrino NLSP which subsequently decays 
to a neutrino and a gravitino. 

An interesting implication of the 3-body kinematics for collider signatures of squarks is a 
reduction in the visible and missing energy of the events, decreasing the efficiency of these 



types of signals to pass standard cuts and making its discovery more challenging. Fig, 14 
contrasts the distributions of missing transverse energy and Ht^ for the two types of decay, 
assuming a squark mass niq = 400 GeV. For the 2-body topology, the neutralino mass 
is chosen to be my,o = 200 GeV, and for the 3-body topology the two invisible final state 
particles have masses m^o = m^o = 100 GeV. 

As a brief illustration of the effects of this different kinematics, we perform a Monte Carlo 



simulation of the two benchmark scenarios introduced in Fig, 14 , and compare their relative 
efficiency to pass a standard jets plus $rp search. Specifically, we consider the baseline 
selection in the multijets plus missing transverse energy search performed by CMS |16] with 
4.98 fb^^. Such selection required events with Ht > 500 GeV, |fr> 200 GeV, and at least 
3 central jets with pt > 50 GeV, among other requirements. Assuming first and second 
generation degenerate squarks, this CMS search constrained the production cross section of 
the 2-body benchmark to be Cqq* < 1.2 pb. This cross section is roughly half of the reference 
QCD cross section at NLO, and therefore the 2-body benchmark is ruled out. The efficiency 
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FIG. 13: Topology of a 3-body squark decay to two invisible sister particles. 



^ Ht is defined as the scalar sum of the px of all jets in the event. 

^ We use MadGraphS for event simulation, Pythia6 for decays, showering and hadronization and PGS4 for 
detector response. 
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FIG. 14: Contrast between the and Ht distributions of 2- and 3-body squark decays. We 
assume rriq = 400 GeV; for the 2-body decay, the mass of the invisible particle is 200 GeV, and for 
the 3-body decay the two invisible particles have mass of 100 GeV. 

for the 3-body topology is roughly 30% of its 2-body counterpart. Naively applying this 
lowered efficiency to the cross section limit, we get a bound aqg* < 4 pb, suggesting that the 
3-body topology is slightly beyond this search's reach. 

Dedicated searches for more challenging SUSY signatures, such as stop and sbottom 
pair production are also impacted by the 3-body decay kinematics. We briefly look into 
one more case, namely, sbottom searches targeting the 2-body topology bi — t- bx^. ATLAS 
performed a sbottom dedicated search [17] with 2 fb~^ of integrated luminosity. The baseline 
event selection required > 130 GeV and two b-tagged jets with > 130 GeV and 
> 50 GeV, among other requirements. Assuming a sbottom mass m^^ = 300 GeV, 
we compared the selection efficiencies for the 2- and 3-body topologies, assuming 50 GeV 
for the mass of the invisible particle in the 2-body decay, and 25 GeV for both invisible 
particles in the 3-body decay. The efficiency of the 3-body topology relative to the 2-body 
topology is ~ 33%. Additional signal regions were defined by placing cuts on the boost- 
corrected contransverse mass |l8], mcT, which have even lower relative efficiencies. The 
ATLAS search excluded this 2-body decay mass point, bounding the sbottom production 
cross section to be cr^^^* < 0.44 pb. Conservatively using the relative efficiency of the baseline 
selection to infer the bound on the 3-body topology, we obtain cr^^j;, < 1.3 pb, which is at 
the borderline of ATLAS' exclusion reach. 

While we have only touched on the phenomenology of such decays, we can see that there 
are clearly interesting results: squarks of 400 GeV and sbottoms of 300 GeV may plausibly 
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be allowed under current LHC limits. The impact of 3-body topologies in LHC searches is 
a rich subject that deserves more careful investigation. We defer this to future work. 

VI. THE COLLIDER PHYSICS OF G-QUARKS 

The phenomenology of the G-quarks is intimately tied into a number of questions: first, 
are we considering a scenario in which there is a preserved sister charge or not? If there 
is a sister charge, is it identified with any of the existing global standard model quantum 
numbers (i.e., baryon or lepton number)? Moreover, if there is a sister charge, and R-parity 
is conserved, is the lightest sister-charged particle (LSiP) also the lightest supersymmetric 
particle (LSP)? 

Let us begin with the simplest case - there is a single vectorlike pair of G-quarks Dg and 
D^g. These fields can interact with the standard model via ^DgD'^ which preserves a U{l)s- 
We could have instead included ^DgD'^ + Ti^DgQ, which gives opposite U{l)s assignments, 
but since Dg falls into a 5 with S^, it is natural to give them the same charge assignments. 

If this sister charge is broken, it can be broken spontaneously or explicitly. If broken 
spontaneously by (p and a vevs, it should be gauged to avoid the presence of a dangerous 
goldstone boson. In this case the production of G-quarks will naturally yield D^^^ — )■ D'^(f){a). 
With (f){cr) — 7- bb we can naturally expect six-jet events with three-jet resonances. CMS has 
performed a search for three-jet resonances in the context of hadronic RPV-gluino decays 
|49j . excluding gluino masses rrig < 460 GeV. The limits from this search on G-quarks 
are much weaker, since the DgDg production cross section is roughly an order of magnitude 
smaller than gluino pair-production. A more quantitative study of the limits on such objects 
can be found in [HO] . 

The G-squarks will mix with SM squarks, and so if they are lighter than the G-quarks, 
they will produce conventional jet+-^2- signals when pair-produced. However, if heavier 
than the G-quarks, and assuming a conventional LSP such as a Bino, they will decay dg — >■ 
Dg + x^ — i- jbbx^. Thus, they will produce six jet + events, with large numbers of b-tags. 
For large G-quark -G-squark splittings, standard jets+^x searches might be sensitive to 
these signatures, although with smaller reach relative to MSSM squarks due to the reduced 
G-squark cross section. (The limits on MSSM squarks usually quoted by LHC searches 
implicitly assume degenerate 1st and 2nd generation squarks, with an effective cross section 
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FIG. 15: (a) Gluino decay topology to dgDg, resulting in final states with 2j + 46 + x^. (b) 
Benchmark spectrum with small splitting between mj and + m^o. 

eight times larger than the cross section of a single squark species). On the other hand, 
for more compressed spectra the missing energy in the events may be significantly reduced, 
rendering high-^j. searches insensitive to such models. However, searches with high-jet 
multiplicity or/and high b-jet multiplicity may be effective for these topologies. Still, the 
cross section may be the limiting factor in sensitivity to G-quarks. 

An interesting scenario to consider is the case of gluino-boosted G-(s) quark production. 
If the G-(s) quarks are lighter than the gluinos with all other squarks heavy, the on-shell 



decay g — t- Dgdg might dominate. This decay topology is illustrated in Fig, 15|^a) - each 
gluino decays to 2j + 46 + x° final states, for a total of 12 jet + events (with large 
numbers of b-tags). Dedicated searches with high b-jet multiplicity could be sensitive to 
such final states, such as the search performed by ATLAS in missing energy and at least 3 
b-jets This search, however, places a strong requirement on missing transverse energy 
{^rp > 160 GeV), which narrows its sensitivity to spectra with large mass splittings between 
G-quarks and G-squarks. In Figjlsj^b) we illustrate a typical spectrum which the ATLAS 
3-bjets-|--^2- search is not sensitive to. In the rest frame of dg, the 3-momentum of its 
neutralino daughter is \p^o\ ~ 62 GeV, and hence this signal has a very low efficiency 
to pass the > 160 GeV requirement of this search. Despite its low missing energy, this 
particular spectrum is distinctive enough to be caught by high jet-multiplicity searches, and 
indeed, through our MC simulations we estimate that this mass spectrum is excluded by 
the ATLAS search in this channel [52]. This search looks for events with up to 9 hard jets, 
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and bounds the gluino cross section in the benchmark of Fig, 15 to be agg < 0.1 pb, which 
is roughly the value of the reference QCD cross section for a 750 GeV gluino. Hence, this 
point is marginally excluded. 

If this sister charge is not broken by and S vevs, there must be a larger gauged sister 
group, I.e., 5f/(2)™ X U{l)f''''\ which leaves some residual U{l)s- In this case, we must 
determine the identity of the LSiP. Assuming it is stable, it is most easily neutral, such as 
the 0, in which case the decay Dg D + LSiP would resemble a standard SUSY squark 
search. Likewise dg — )■ DgX^ D + + LSiP could be discoverable as well with enough 
luminosity. 

The sister charge can be broken explicitly if the sister group is 5f/(2)™xf/(l)|'°''''\ I.e., 
we can write ?7(1)e breaking terms (pDgD'^ + (j)DgD^ + "EdDgQ. In this case, both G-quarks 
decay is as in the broken U{1) case above, with multijets and many b-tags. 

Another interesting case is when the U{l)s is identified with some SM global charge. As we 
discuss in the appendix, the operator UD'^D'^ naturally identifies U{l)s with baryon number, 
but does not break R-parity. In this case, while the G-quarks without sister charge would 
decay into three jets, the sister-charged G-quark would decay — uDs U + D + bb + XQ-i 
yielding eight-jet plus $rp signals, and sister charged G-squarks would decay dg*^ — t- uD'g — >• 
udbh, i.e., we would have eight jet events, with four jet subresonances. (Here, the G-squarks 
of different sister charge have different R-parity under the intact R-parity after 0, S acquire 
vevs.) 

In summary, the presence of G-quarks can radically alter our expectations for phe- 
nomenology. In particular, if the G-quarks are lighter than the gluinos, it can dramatically 
enhance the number of jets (in particular b-jets) in the event, and suppress the signal. 



VII. CONCLUSIONS 



If there is new physics at the weak scale, the Higgs boson and its properties are the 
first indication of what that may be. The value of 125 GeV constrains a wide range of 
supersymmetric models, and may be pointing us to a new contribution to the Higgs quartic. 
Models such as the NMSSM can raise it, but at the cost of including pure singlets, and only 
with very large couplings. 

In the context of the NMSSM, however, it is possible to replace Hd with a new field S^, 
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a "Sister Higgs," and still yield the usual tree-level corrections to the Higgs mass. Such a 
field need not couple to ordinary matter, and, indeed, would not if it carries a charge under 
some new sister group Gs- If that group is gauged, new contributions to the NMSSM-like 
superpotential term can drive it large, easily boosting the Higgs mass up to 125 

GeV, or beyond. 

The presence of this new Higgs has a number of consequences: as with a Type I 2HDM, 
a light charged Higgs is now possible. The overall consequences allow a modification to the 
Higgs properties, such as a boosted VBF 77 signal, as well as a boosted inclusive 77 signal, 
with suppressed decays to fermions [53J. 

Once we have expanded the theory to include the sister Higgs, we could also expect yet 
more Higgses when Gs is non-Abelian. With these fields we can write down the R-parity 
violating operator Hd^dE'^ once Gs is broken. Because this operator couples only to right 
handed fermions, it is far less constrained that traditional lepton-number violating RPV 
operators. While it changes the phenomenology of SUSY decays, if the finals state contains 
r's, the current search strategies should find it. At the same time, it motivates searches for 
leptoquarks, with squarks decaying to but importantly without an associated missing 
energy signal (as there is no accompanying channel with a neutrino). 

The Z' in these models is natural, and can be quite light without significant constraints 
from the p parameter or searches for dilepton resonances. Corrections to Zh production can 
be sizeable at the Tevatron but are not expected to be large at the LHC without additional 
couplings to SM fermions (although such corrections could come from mixing with the G- 
quarks at the cost of fiavor constraints). 

Appeals to grand-unification motivate us to consider the presence of colored fields which 
are charged under Gs as well. Unlike conventional fourth-generations, these vectorlike quarks 
decay dominantly through emission {or (p — h mixing induced h emission) . These G-quarks 
can have important effects on phenomenology. In particular, if they are lighter than the 
gluino, they can open up new decay channels with lower missing energy and more jets than 
conventional gluino decays. The phenomenology of G-quarks (and G-squarks) is very rich 
and we have only begun the discussion in this work. 

In the presence of a sister charge, the lightest sister-charged particle (LSiP) would be 
stable. If it is a Dirac fermion, constraints from direct detection searches would have excluded 
it in all but the narrowest corners of parameter space where the WIMP is very under 
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abundant or is very light. It may thus be most natural that symmetry should be broken, at 
least to a Z2, or become identified with B or L so that the LSiP can decay into it. 

Regardless, if the LSP also carries sister charge, squarks may only decay via three-body 
decays. Such decays have reduced missing energy and Ht, and are consequently less con- 
strained by existing iets+$x searches. 

In summary - very complex phenomenology can arise by the simple extension of the MSSM 
by the presence of G^-charged sister-Higgs fields and their related G-quarks. A low energy 
theory that is not the MSSM, nor even the NMSSM, but one expanded with additional Higgs 
fields easily provides for a Higgs mass as large as 125 GeV and a profoundly changed LHC 
phenomenology. As more data accumulate, we shall see soon if such a rich Higgs sector is 
realized in nature. 
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